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ABSTRACT 

Aims. We discuss the properties of the progenitors of core collapse supernovae of type lb and Ic and of long soft gamma 
ray bursts, as they can be deduced from rotating stellar models of single stars at various metallicities. 
Methods. The type of the supernova progenitor was determined from the surface abundances at the pre-supernova 
stage. The type of the supernova event was obtained from the masses of hydrogen and helium ejected at the time of 
the core-collapse supernova event. 

Results. We find that the minimum amount of helium ejected by a core-collapse supernova (of whatever type) is around 
0.3 Mq. There is no diff'erence between the WC and WO stars in the ejected masses of helium, CNO elements, and 
heavy elements. Also no difference is expected between the chemical composition of a WC star resulting from a normal 
or a homogeneous evolution. The progenitors of type lb supernovae are WNL, WNE, or less massive WC stars. Those 
of type Ic are WC and WO stars. WO stars are produced in a limited mass range (around 60 Mq) and only at low 
metallicity (for Z <~ 0.010) as already found. The WO stars are the progenitors of only a small fraction of type Ic. 
Present stellar models indicate that, at solar metallicity, there is about 1 type lb supernova for 1 type Ic, and this ratio 
rises to 3 type Ic for 1 type lb SN at twice solar metallicity. At this metallicity, type Ic's are more frequent than type 
Ib's because most massive stars that go through a WNE stage evolve further into a WC/WO phase. Current models 
can account for the observed number ratios SN Ib/SN II and SN Ic/SN II and for their observed variation with the 
metallicity. In case no supernova occurs when a black hole is formed, single-star models can still account for more than 
half of the observed (SN Ib-fSN Ic)/SN II ratio for Z > Zq. For the gamma ray burst rate, our models produce too 
large a number for such an event, even if we restrict the progenitor to the WO stars. This confirms that only a fraction 
of the WC / WO stars evolve toward gamma ray burst event, most likely those arising from stars that were initially 
very rapid rotators. 

Key words, stars: Wolf-Rayet, supernovae: general, gamma rays: bursts 
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On 
O 



„ , novae. Such a trend is reproduced well by rotating single- 

Type lb supernovae are core-collapse supernovae whose ^^^^ m odels accounting fo r metallicity-dependent stellar 

spectrum shows no hydrogen Imes. The spe ctra of type Ic ^-^^^ (|Mevnet fc Maedeil [200l . On the other hand, a 

s^ no hydrogen an dMrum hues (see e.g. [ Wheeler et al. | ^^jj ^^^^^ ^^^^^^.^ explaining this tr end comes from close 

lim i Nomoto et al. | 11994). The progenitors of these core binar y-star evolution (se e for examplelPodsiadlowski^ialJ 

. collapse supernovae are thus believed to be stars stripped of ^IvanbevereT^TO lEldridee et al. 2008.) . In this 

their original H-rich envelope for type lb sand also of their ^^^^^^^^^ ^^e hydrogen-rich envelope is removed either 

He-rich envelope for type Ic s Progenitors are therefore ^^^^^ ^ ^^^^^ ^^^^ overflow process or during a com- 

Wolf-Rayet stars. The observed frequency at solar metal- ^^^^^ ^ producing a WR star, 
licity 01 type Ibc supernova e is about 20yo the fre quency of 

type II supernovae (see e.g. ICappellaro et Zl fT999l which ^^^^ preseiit work, we takejie opportunity of recent 

represents a significant fraction of all core collapse super- observations by| Prieto et al. | (|2008D from which separated 

novae. In four cases, the typical spectrum of a type Ic super- frequencies for type lb and Ic can be deduced to get a step 

nova has been observed associated with a long so ft gamma ^^''t^^er. The questions we want to address are: how many 

ray burst (GRB) event (|Wooslev fc BloomllS ). indicat- t^P^ and Ic supernovae normalized to the number of 

ing a privileged link between type Ic's and the most pow- ^^P^ " supernovae can we expect from single-star mod- 

erful supernova explosions observed in the Universe. That ^Is at different metallicities, how the numbers derived from 

type Ibc represents a significant proportion of core collapse smgle-star models compare with the observed values, and 

supernovae and the link between type Ic's and the long soft ^^^^^^ change no supernova event occurs if ev- 

GRBs justify the study of the evolution leading to such ^ry time a black hole forms. Answers to those questions 

events ^^^^ provide some hints as to the importance of the single 

In this context, a very interesting feature observed by ^^^^ scenario in explaining the evolution leading to the type 

iPrantzos fc Bo"iM (IM) is the increase with the metal- ^^c supernovae events and will provide useful quantities for 

more comparisons with similar outputs obtained from close 

Send offprint requests to: Cyril Georgy binary evolution scenarios. 
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Other questions will be addressed, for instance, the ab- 
sence of H and He lines does not necessarily imply that the 
associated elements are completely absent from the ejecta. 
It only tells us that the physical conditions (the chemi- 
cal composition being only one of them) are such that no 
lines of these elements are formed at the time of the super- 
nova event. Actually, the quantities of hydrogen and helium 
ejected by these two types of supernovae are not known. In 
this paper, guided by stellar models for single stars, we es- 
timate the minimum quantities of hydrogen and of helium, 
we may expect in type lb and type Ic supernovae. We also 
indicate the relations of filiation between type lb, Ic super- 
novae and the WNL, WNE, WC, and WO stars. 

The second section of this paper briefly recalls the main 
physical ingredients of the stellar models used in this study. 
Based on these models, we discuss in Sect. |3]the nature of 
the supernova progenitors, and in Sect.|3]the chemical com- 
position of the supernova ejecta. In Sect.[5l we discuss the 
link between the chemical composition of the ejecta and the 
nature of the supernova. The nature of the stellar remnant 
is discussed in Sect. |6l Theoretical rates for different su- 
pernova types are estimated in Sect. [7] and compared with 
observations. Section [5] discusses the implications for the 
long soft GRB progenitors and finally conclusions are pre- 
sented in Sect. [H 

2. Stellar models 

All the rotating s t ellar model s we use here come from 
iMevnet k Maedeil (|200a Hool) (hereafter papers X and 
XI). A complete description of the physics used can be 
found in these two papers. We just recall here some im- 
portant characteristics. 

— All the models considered in this work are computed 
with an initial equatorial velocity of 300kms~^. This 
initial velocity implies time-averaged equatorial veloci- 
ties between 160 and 250 km depending on the ini- 
tial mass and metallicity (see tables 1 in papers X and 
XI). These values are in the range of observe d values 
for OB stars; for instance. iHuang fc GiesI ()2006D present 
projected rotational velocities for 496 OB stars belong- 
ing to 19 young galactic clusters with estimated ages be- 
tween 6 and 73 Myr. Mean usini values of 139, 154, and 
151 km s~^ were obtained for groups of 09.5-B1.5, B1.5- 
B5.0, and B5.0-B9.0 type stars. These authors derived 
the underlying probability distribution for the equato- 
rial velocities v and obtained a peak at 200 km s~^. 
iDufton et al.l ()2006D obtain a peak of v at 250 km s^^ 
with a full width half maximum of approximately 180 
km s~^ for the unevolved targets in the galactic clusters 
NGC 3293 and 4755. 

— Mass loss rates are those of IVink et all (|2000l l200l| ). 
w hich take wind bistability into account. The rates 
of Ide Jager et al.l (|l988f ) are used outsid e the domai n 
of apphcation of the rates of IVink et aH (|2000l l200lh . 
During the WR phase, the mass loss rates are those of 
iNugis fc Lamersi (j200a ). These mass loss rates, which 
account for the clumping effects in the winds, are lower 
by a factor of 2-3 than t he ones used in our previous 
non-rotating stellar grids (jMevnet et al.l[T993 ). 

— The dependence of mass loss on rotation 
(|Maeder fc Mevnetl 120001 ) was taken into account, 
as is the dependence on the metallicity. During the 



non-WR phases of the present models, we assumed that 
the mass loss rates depend on the initial metallicity 
as M(Z) ^ (Z/Z p,V/^M(Zp,) (|Kudritzki fc Puk.2000l : 
IVink et all 1200 lH . 

Models were computed for four metallicities: Z — 0.004, 
Z = 0.008, Z = 0.020 (standard) and Z = 0.040. The initial 
masses considered are indicated in Table [Tj All the models 
were computed up to the end of the core helium-burning 
phase. The very short durations of the advanced phase im- 
ply that the star at the end of the core He-burning phase 
has already reached its final mass. Moreover, the decou- 
pling between the rapid evolution of the core governed by 
the neutrino emission and the slow evolution of the enve- 
lope does not allow the envelope to change a lot during the 
advanced phases. Thus we consider that the properties of 
the outer layers of our models obtained at the end of the 
core He-burning phase are very near the one that would 
have obtained had we pursued the computation until the 
presupernova stag43- 

In previous papers ( s ee e.g. Mevnet fc Maedeil 120001 : 
iMaeder fc Mevnetl 1200 ll : iMaeder et all l2008l Paper X; 
Paper XI), we compared the outputs of the rotating models 
to many observed features of massive stars and showed that 
the models accounting for the effects of rotation gave much 
better agreement for fitting observed features such as 

— surface enrichments, 

— the blue to red supergiant ratio in the SMC, 

— the variation with metallicity in the number ratio of WR 
to 0-type stars, 

— the observed number ratio of WN to WC stars in the 
LMC and SMC, 

— the variations with metallicity in the number ratio of 
type Ibc to type II supernovae, 

— the existence of Wolf-Rayet stars showing both at their 
surface H- and He-burning products. 

In papers X and XI we discussed the last four points 
above in details. Here, as indicated in the introduction, we 
focus on the nature of the progenitors for the type lb and 
type Ic SNe, taken separately. 

3. Nature of supernova progenitors 

The nature of the supernova progenitors was determined 
by the surface composition of the star at the presupernova 
stage. The following classification scheme was adopted: 

— All the stars not ending their lifetime as a Wolf-Rayet 
star end their nuclear life as a blue or a red supergiant. 
They are hereafter called supergiant (SG). 

— The star is considered to be a WR star when log T^g > 
4.0 and the mass fraction of hydroge n at the surface Xfj 
is ijnferior to 0.4 (same criterion as in lMevnet fc Maedeil 



— Wolf-Rayet stars with hydrogen at their surface {Xs > 
10-^) are WNL stars. 

— Wolf-Rayet stars without hydrogen on their surface 
{Xs < 10^^) and with surface carbon abundance in- 
ferior to nitrogen abundance are WNE stars. 



We checked using the solar metallicity models of 
iHirschi et all (120041 ) . which were pursued until the presupernova 
stage, that indeed this is the case. 
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~ Wolf-Rayet stars with carbon abundance on the sur- 
face superior to nitrogen abundance are WC or WO 
stars. To distinguish b etween WC and WO, we follow 
ISmith fc Maederl (jl99lh : if the ratio of (in number) 
is less than 1, we have a WC, otherwise a WO. 

A star enters the WR phase as a WNL, then may evolve 
through the sequence WNL^ WNE ^ WC WO. A star 
can end its evolution at any of these stages. We do not 
consider the WN/WC phase here, which can occur as a very 
short transition phase between the WNE and WC phases. 

Table [1] gives for the 4 metallicities considered here the 
initial mass of the star Mini, the mass at the end of the 
helium-burning phase Mend He, the mass of the He-core, 
Mho, defined as the mass interior to which the mass frac- 
tion of helium (or of its products after fusion) is higher than 
0.75, the mass of the carbon-oxygen core, Mcoj defined as 
the mass interior to which the mass fraction of carbon plus 
oxygen is higher than 0.75, th e mass of the remnant A/icm 
(given bv lHirschi et al.|[2005f) . the type of the progenitor, 
the SN type produced (see Sect.[S]), and the remnant type 
(see Sect. 

In Fig.[l] we show how the initial mass ranges of various 
supernova progenitors vary as a function of the metallicity. 
To build this figure we had to adopt a method to determine 
the mass limits between the different zones. To illustrate, 
let us consider one specific case. At Z = 0.02, from Table[T] 
we know that the 25 Mq stellar model ends its lifetime as 
a WNL star and the 40 Mq as a WC star, but we know 
neither at which mass the transition from WNL to WC 
star occurs nor whether some stars will end their life in the 
WNE stage. To find the transition between the WNL and 
WNE stage, we deduce from the chemical structure of our 
25 Mq at the end of its evolution, the mass that should be 
removed in order for the star to become a WNE star (in 
that case 0.35 M0 or 1.4% of its total initial mass). From 
our 40 Mq stellar model at the end of its evolution, we 
deduce the mass it has lost since its entry into the WNE 
phase (1.53 Mq or 6.1% of its initial mass). Assigning, by 
convention, a minus sign to this last quantity (mass lost) 
and a positive sign to the former one (mass which has to 
be lost), plotting these two quantities in a plane "mass- 
to be/or- lost" versus "initial mass", the abscissa of the 
intersection of the line connecting these two points with 
the horizontal line at ordinate gives the initial mass for 
which the cndpoint of the evolution just corresponds to the 
stage where "no mass" has to be lost to enter into the WNE 
phase and where no mass has been lost since the entry into 
that phase. Thus this is the minimum initial mass to end 
its evolution as a WNE star. Using similar procedures for 
all the limits, we obtain the regi ons shown in F i g. [1] N ote 
also that we used the results of lEkstrdm et al.] (120081 ) for 
complementing the figure with Pop III stellar models. 

Looking at that figure, we can make the following re- 
marks. 

— As expected, the lower mass limit for having a WR pro- 
genitor decreases when the metallicity increases. We see 
that the dependence on the metallicity of this limit is 
much stronger at low Z than at high Z. 

— The mass range of stars ending their life as a WN star 
is relatively narrow. Below this mass range, the stars do 
not succeed in entering the WR phase; above it, mass 
loss rates are efficient enough for allowing the star to 
evolve further into the WC or WO phase. 



— We note that at metallicities higher than about 0.02, 
the WN progenitors are more or less equally distributed 
among WNL and WNE stars, while for metallicities be- 
low 0.01 all WN progenitors are WNL stars. Actually 
present stellar models predicts such a narrow range of 
masses for WNE progenitors below Z = 0.008 that 
we neglected it. This can be explained in the following 
way: first let us point out that WNE stars can be pro- 
duced only when the star has evolved beyond the core 
H-burning phase. Indeed, regions without hydrogen are 
present in the star only after the main-sequence phase 
has ended. At high metallicity, because of strong stellar 
winds, the outer part of the H-burning core is uncov- 
ered at an early stage when the size of the He-burning 
core is still small and has not yet transformed a signifi- 
cant part of the former H-burning core into carbon and 
oxygen. This allows the intermediate region between the 
He-burning core and the H-burning shell to have a large 
extension in mass. Thus, the WNE phase has a signi f- 
icant duration (see Fig. 9 of iMevnet k Maedeil 12005V 
which allows a narrow range of initial masses for which 
the endpoint of the evolution occurs when the star is 
in the WNE phase. At low metallicity, the outer part 
of the H-burning core is uncovered at later evolutionary 
stages when the He-burning core has grown in mass, 
reducing the extent in mass of the region with chemi- 
cal composition typical of the WNE phase. This reduces 
the duration of the WNE phase and also the initial mass 
range of stars ending their evolution in that stage. 

— Above 30 Mq and for metallicities higher than 0.008, all 
stars end their life as WC/WO stars. For metallicities 
below 0.008, the minimum initial mass for stars ending 
their lifetime as WC stars rapidly increases when the 
metallicity becomes lower. 

— The WO progenitors region appears as an "island" in 
the "WC" sea, located at low me tallicity. As first ex- 
plained bv lSmith fc Maede^ ()1991[ ). the low metallicity 
position of this island comes from the fact that, in order 
to form a WO star, the He-burning core must be uncov- 
ered at a very advanced stage of its evolution when a lot 
of helium has been converted into carbon and oxygen. 
This is more easily realized at low metallicity where the 
mass loss rates are lower, hence the core is revealed at 
a later stage. We also note that the WO progenitors do 
arise neither from the most massive nor from the less 
massive stars having their He-burning core uncovered. 
This reflects the delicate interplay between two counter- 
acting effects. On one hand, in order for a star to become 
a WO star, strong enough stellar winds are necessary; 
on the other hand, as explained above, the stellar winds 
should not be too strong, otherwise no WO stars but 
instead WC stars are formed! 



4. Chemical composition of the ejecta 

The masses of hydrogen and helium ejected at the time of 
the explosive event are the key quantities for determining 
the type of the supernova event. These quantities are given 
for each of our stellar models in Tabled which indicates for 
each metallicity and initial mass considered the total mass 
ejected at the time of the supernova explosion (column 2), 
the mass of hydrogen, of helium, of carbon, nitrogen and 
oxygen, and of heavy elements (columns 3 to 8). Heavy 
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Table 1. For all our models, various quantities are given here (see text for details). 



Mini 


Mend He 




Mco 




Prog. 


SN 


Rem. 


Mini 


Mend He 


Mne 


Mco 


M-en, 


Prog. 


SN 


Rem. 








Z = 0.004 














Z = 0.020 








12 


11.8 


3.74 


1.78 


1.4 


SG 


II 


NS 


12 


10.5 


3.81 


1.94 


1.4 


SG 


II 


NS 


15 


14.1 


5.01 


2.84 


1.7 


SG 


II 


NS 


15 


10.2 


5.97 


3.61 


1.7 


SG 


II 


NS 


20 


18.0 


7.45 


4.79 


2.0 


SG 


II 


NS 


20 


11.8 


9.01 


6.62 


2.0 


SG 


II 


NS 


25 


20.0 


9.95 


7.06 


2.5 


SG 


II 


NS 


25 


11.3 


11.33 


8.98 


2.5 


WNL 


II 


NS 


30 


18.9 


14.16 


11.14 


2.9 


SG 


II 


BH 


40 


12.7 


12.70 


11.75 


3.0 


WC 


Ic 


BH 


40 


22.3 


21.59 


17.16 


3.5 


WNL 


II 


BH 


60 


14.6 


14.60 


13.67 


3.6 


WC 


Ic 


BH 


60 


28.5 


28.50 


28.46 


4.4 


WO 


Ic 


BH 


85 


12.3 


12.30 


11.22 


3.0 


WC 


Ic 


BH 


120 


17.2 


17.20 


17.18 


3.4 


WC 


Ic 


BH 


120 


11.3 


11.30 


10.40 


2.8 


WC 


Ic 


BH 








Z = 0.008 














Z = 0.040 








20^ 


15.5 




4.85 


2.1 


SG 


II 


NS 


20 


9.2 


9.05 


6.84 


1.4 


SG 


II 


NS 


25^ 


14.0 




6.92 


2.7 


SG 


II 


NS 


25 


9.6 


9.60 


8.19 


1.7 


WC 


lb 


NS 


30 


12.1 


12.10 


11.10 


2.9 


WC 


Ic 


BH 


40 


11.4 


11.40 


10.48 


2.8 


WC 


Ic 


BH 


40 


17.4 


17.40 


17.23 


3.5 


WC 


Ic 


BH 


85 


7.2 


7.20 


6.07 


2.2 


WC 


Ic 


NS 


60 


16.4 


16.40 


16.44 


3.3 


WO 


Ic 


BH 


120 


7.1 


7.10 


6.10 


2.2 


WC 


Ic 


NS 


120 


13.4 


13.40 


13.24 


3.0 


WC 


Ic 


BH 



















Homogenous model at Z — 0.002 

60^ 27.8 - - 4.0 WC lb BH | 

(1) Non-rotating models from Mevnet et al. (1994) used to estimate the minimum initial mass of stars ending their lifetime 
in the WNL phase at that metallicity. These models were computed with artificially enhanced mass loss rates, which 
somewhat compensates for the neglect of the effects of rotation. 

(2) Model computed up to the end of silicon burning. 

Note: All masses are given in [Mq]. The SN type is given for mg'" = 0, m^"" = 0.6 and su pposing that a SN occurs even when 
a BH is formed (see text). The remnant masses used are those given bv lHirschi et al.l (|2005l ). except for the homogenous models 
(see text). 



elements are the sum of the mass ejected under the form of 
elements heavier than hydrogen an d helium. For ca rbon and 
oxygen, we used relations given bv lMaedeil (1X9921 ) between 
the mass of the CO core and the ejected mass of these two 
elements at the time of the supernova event. 

From the results shown in Table [H we can deduce the 
following trends. 

— Most of supergiant progenitors eject more than 70% of 
their ejecta in the form of hydrogen and helium. At any 
given metallicity, the proportion of the ejecta under the 
form of hydrogen and helium is greater in lower initial 
mass stars. 

— The 20 Mq stellar model aX Z = 0.04 is a special case 
of supergiant progenitor in the sense that it only ejects 
small amounts of hydrogen and helium. This model is 
a blue supergiant (log Ted = 4.338) at the end of its 
stellar lifetime and is actually on the verge of becoming 
a Wolf-Rayet star. It should lose only 0.09 Mq to un- 
cover layers with a hydrogen mass fraction inferior to 
0.4. Probably a slightly more massive star would show 
an evolutionary path of the type 0-red SG-blue SG- 
WR. Let us mention here that the single-star evolu- 
tionary scenario predicts a strong anti-correlation in the 
numb er of RSG and o f WR stars in single-aged popula- 
tions (|Leitheredll999| ). Close binary evolution may pro- 
duce WR stars more easily from lower initial mass stars 
(jPodsiadlowski et al.lll993 ) and thus explain the simul- 
taneous presence of RSG and WR stars even in stel- 
lar populations with age spread inferior to a few Myr. 
It would be interesting to obtain from well-observed, 



single-aged stellar populations quantitative constraints 
on the occurrence of RSG and WR stars and to compare 
them with population synthesis models accounting for 
single and close binary evolution. 

— Our computed WNL progenitors eject less than 1 Mq 
of hydrogen, but more than 1 Mq of helium (more gen- 
erally every time some hydrogen is present in the ejecta, 
the mass of ejected helium is above ~2 Mq). They eject 
greater quantities of carbon and oxygen than supergiant 
progenitors. About 3/4 of the mass of the ejecta is in 
the form of heavy elements. 

— WC and WO progenitors eject more than 90% of their 
ejecta in the form of heavy elements. In these cases, the 
mass of ejected helium is between 0.31 and 0.54 Mq. 

— We see that the minimum helium content in the ejecta of 
a core collapse supernova is about 0.3 Mq. This means 
in particular that type Ic SNe should at least contain 
this minimum amount of helium in their ejecta. This 
agrees with the resul ts of n on rotating models presented 
in E ldridge fc Tout) (j2004f ). and with binary models of 
Eldridge (private communication). There is to date no 
clear measurement of the mass of heli um contained in 
the ejecta of type Ic SNe (although see lElmhamdi et al.l 
2006.) . It would be very interesting to find a correlation 
or an anticorrelation between the mass of ejected he- 
lium and the mass of another element produced in the 
same region of the star that would be more easily de- 
tectable. Typically, such correlations exist between the 
ejected mass of He and of N, as can be seen in Fig. [31 
We see that the correlation depends on the metallic- 
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Table 2. Initial mass Mini, total released mass M, 



O (Mo) and of heavy elements {M\ 



heavy J 



1, mass and fraction of ejected H (A/h), He (A/hc) 
with all masses in Mq. 



C (A^c), N (A/n), 



Mini 



mc 



rriN 



mo 



mil. 



0.004 



12 


10.4 


5.52 (53%) 


4.30 (41%) 


0.15(1.4%) 


1.7 


10-2 (0.16%) 
10-2 (0.12%) 
10-2 (0.11%) 
10-2 (0.12%) 
10-2 (0.13%) 
10-2 (0.07%) 


8.9 ■ 10-2 (0.86%) 


0.58 (5.6%) 


15 


12.4 


6.11(49%) 


4.94 (40%) 


0.25 (2.0%) 


1.5 


0.42 (3.4%) 


1.35(11%) 


20 


16.0 


6.57(41%) 


6.35 (40%) 


0.49 (3.1%) 


1.8 


1.41 (8.8%) 


3.08 (19%) 


25 


17.5 


5.48 (31%) 


7.09 (41%) 


0.88 (5.0%) 


2.1 


2.64(15%) 


4.93 (28%) 


30 


16.0 


1.58 (9.9%) 


5.75 (36%) 


1.45 (9.0%) 


2.1 


5.07(32%) 


8.67 (54%) 


40 


18.8 


0.45 (2.4%) 


4.35 (23%) 


2.11(11%) 


1.4 


9.93 (53%) 


14.0 (75%) 


60 


24.1 


0.00 (0.0%) 


0.41 (1.7%) 


2.16(9.0%) 


1.4 


■ 10-" (0.0%) 


14.44 (60%) 


23.69 (98%) 


120 


13.8 


0.00 (0.0%) 


0.48 (3.5%) 


1.73(13%) 


0.00 (0.0%) 


6.92 (50%) 


13.32 (97%) 



0.008 



20^ 


13.4 


9.70 (73%) 


3.55 (27%) 


2.4- 10-2 (0.18%) 
1.6 ■ 10-2 (0.14%) 


7.5 ■ 10-^ (0.06%) 
1.5 ■ 10-2 (0.14%) 
1.8 • 10-=^ (0.02%) 
3.0 ■ 10-^ (0.0%) 


5.6 • 10-2 (0.42%) 
4.3 • 10-2 (0.38%) 


0.11(0.8%) 


251 


11.3 


7.83 (69%) 


3.40 (30%) 


0.91 (0.8%) 


30 


9.2 


0.00 (0.0%) 


0.54 (5.9%) 


1.06(12%) 


3.30 (36%) 


8.66 (94%) 


40 


13.9 


0.00 (0.0%) 


0.31 (2.2%) 


1.75(13%) 


6.98 (50%) 


13.59 (98%) 


60 


13.1 


0.00 (0.0%) 


0.49 (3.7%) 


1.68(13%) 


0.00 (0.0%) 


6.48 (49%) 


12.61 (96%) 


120 


10.4 


0.00 (0.0%) 


0.44 (4.2%) 


1.36(13%) 


0.00 (0.0%) 


4.55 (44%) 


9.96 (96%) 



12 


9.1 


4.13 (45%) 


4.20 (46%) 


15 


8.5 


2.53 (30%) 


3.83 (45%) 


20 


9.8 


1.40 (14%) 
3.78 ■ 10-' (0.0%) 


3.51 (36%) 


25 


8.8 


2.04 (23%) 


40 


9.7 


0.00 (0.0%) 


0.53 (5.5%) 


60 


11.0 


0.00 (0.0%) 


0.54 (4.9%) 


85 


9.3 


0.00 (0.0%) 


0.53 (5.7%) 


120 


8.5 


0.00 (0.0%) 


0.49 (5.8%) 



Z = 0.020 

0.18(2.0%) 
0.32 (3.7%) 
0.74 (7.6%) 
1.13(13%) 
1.21(12%) 
1.43 (13%) 
1.11(12%) 
0.96(11%) 



5.4- 
5.0- 
4.7- 
2.6- 
3.8 ■ 
4.4 ■ 
3.7- 
3.4- 



10-2 (0.59%) 
10-2 (0.59%) 
10-2 (0.48%) 
10-2 (0.29%) 
10-" (0.0%) 
10-" (0.0%) 
10-" (0.0%) 
10-" (0.0%) 



0.18(1.9%) 
0.84 (9.9%) 
2.22 (23%) 
3.50 (40%) 
3.77(39%) 
4.86 (44%) 
3.43 (37%) 
2.94 (35%) 



0.77(8.5%) 
2.14 (25%) 
4.89 (50%) 

6.76 (77%) 
9.17(95%) 
10.46 (95%) 

8.77 (94%) 
8.01 (94%) 



Z = 0.040 



20 


7.8 


7.1 ■ 10-2 (0.91%) 


1.96 (25%) 


0.77 (9.9%) 


4.62 ■ 10-2 (0.59%) 


2.29 (29%) 


5.77 (74%) 


25 


7.9 


0.00 (0.0%) 


0.51 (6.5%) 


0.65 (8.2%) 


0.00 (0.0%) 


1.88 (24%) 


7.39 (94%) 


40 


8.6 


0.00 (0.0%) 


0.41 (8.3%) 


1.00(12%) 


0.00 (0.0%) 


3.07 (36%) 


7.89 (92%) 


85 


5.0 


0.00 (0.0%) 


0.48 (9.6%) 


0.33 (6.6%) 


0.00 (0.0%) 


0.89 (18%) 


4.52 (90%) 


120 


4.9 


0.00 (0.0%) 


0.44 (9.0%) 


0.33 (6.8%) 


0.00 (0.0%) 


0.90 (18%) 


4.46 (91%) 








Homogenous model at Z = 0.002 






602 


23.8 


0.00 (0.0%) 


0.64 (2.7%) 


2.33 (9.8%) 


3.66 • 10-^ (0.0%) 


15.56 (65%) 


23.12 (97%) 



(2) Model computed up to the end of silicon burning. 



Table 3. Mass limits between type II and type lb SNe {M^J^^), 
(see text). 



and between type lb SNe and Ic SNe {mI^J^^'') in Mq 



z 


mass limit between type II 
mg'" = 0.00 mg'" = 0.25 


- type lb SNe 
mg'" = 0.50 


mass 


limit between type lb 
= 0.4 mSi" = 0.5 


- type Ic SNe 
= 0.6 


0.004 


54.0 


48.1 


39.4 






59.4 


58.6 


0.008 


30.0 


29.9 


29.7 


35.7- 


- 49.1 


31.6 


30.0 


0.020 


25.1 


24.2 


23.1 






110.2 


39.1 


0.040 


20.5 


20.5 


20.5 






26.2 


24.7 



ity so that to use such theoretical guideline for the He- 
content, some measurement of the metallicity is needed 
(such correlations would allow an idea of the metallicity 
in case the amounts of helium and of nitrogen in the 
ejecta can be measured!). To advance on the question 
of the He-content in type Ic supernova ejecta, we would 
need similar correlations but of course with another el- 
ement, since in type Ic's the amount of nitrogen is zero. 
Probably ^^Ne would be an interesting candidate, since 
it traces the He-rich regions of the star. The difficulty 



would then be to distinguish the ^^Ne abundances from 
that of ■^'^Ne in the spectrum analysis. In a future work, 
we shall study this question. 
— There is no difference in the chemical composition of 
the ejecta of the WC and the WO stars. 

5. The supernova types 

As already emphasized in the introduction, the conditions 
for producing type lb and type Ic supernovae are based 
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Table 4. Mass ranges for various supernova progenitors, supernova types, and compact remnants at different metallicities 
(in Mo). 



Z SG WNL WNE WC WO Type II Type lb Type Ic NS BH 



Reference case 

0.004 8 - 32 32 - 54 - 54 - 59; 114 - 120 59 - 114 8 - 54 54- 59 59 - 120 8 - 29 29 - 120 

0.008 8 - 25 25 - 30 - 30 - 41; 113 - 120 41 - 113 8 - 30 - 30 - 120 8 - 29 29 - 120 

0.020 8 - 23 23 - 26 26 - 29 29- 120 - 8 - 25 25 - 39 39 - 120 8 - 35 35 - 120 

0.040 8 - 20 20 - 20.5 20.5 - 22 22 - 120 8 - 20.5 20.5 - 25 25 - 120 8 - 120 

(1) At this metallicity, only the 40 Mq model produces a BH. 

Note: The SN type is given for ttih'" = 0, rriHo" ~ 0.6 and supposing that an SN occurs even when a BH is formed (see text for 

more details). 




8 10 20 40 60 80 100 



Mini / Ms 

Fig. 1. Ranges of masses of different types of supernova 
progenitors at different metallicities. For each area, the type 
of progenitor is indicated in the figure (see text). The shad- 
ing on the right indicates the area where formation of a 
black hole (BH) is expected ; elsewhere, the remnant is a 
neutron star (NS) (color figure available onhne). 



on the presence/absence of H, He lines in the supernova 
spectrum. Here we adopt a criterion based on the quanti- 
ties of hydrogen and of helium ejected at the time of the 
supernova event to link a given chemical structure of the 
supernova progenitor to a given supernova type. Let us call 
m^™ the minimum quantity of hydrogen that should be 
present in the ejecta for the spectrum of the supernova to 
present H lines and , the minimum quantity of helium 
that should be present in the ejecta for the spectrum of 
the supernova to present He lines. In that case a given ini- 
tial mass star M presenting at the end of its evolution a 
quantity mn^M) of hydrogen in its envelope, respectively 
7TiHc(-^^)j of helium will produce 

— a type II supernova if rnH(Af) > mg'"; 

— a type lb supernova if m}i{M) < m^™ and r7iHc(Af) > 




8 10 20 40 60 80 100 



Mini / M,, 

Fig. 2. Same as Fig. [T] for various types of SNe. The SN 
type is given for to^'" — 0, mg^^" = 0.6 and supposing that 
an SN occurs even when a BH is formed (see text). 

— a type Ic supernova if mH{M) < mg'" and mHc(Af) < 

min 
'"He • 

The quantities m^™ and to^^" are not known a priori. 
In the case of m^'" = 0, we obtain the limiting masses, 
M^^^^^, indicated in the second column of Table 4, where 
stars with initial masses below M^^^^ explode as type II 

supernovae and stars with masses above M^^^^ explode as 
type lb or type Ic supernovae. We see that, when the metal- 
licity becomes higher, this mass limit decreases as expected. 
Interestingly also, we see that adopting m^™ — 0.25, or 0.5 

Mq, pushes m|j'j^^'^ towards lower values, thereby decreas- 
ing the mass range for type II supernovae and increasing 
it for type Ibc supernovae. Therefore m^'" = leads to a 
lower limit for the number fraction of type Ibc to type II su- 
pernovae. Except at the metallicity Z = 0.004, the changes 
remain quite modest. In the following we adopt the value 
m^'" = as our reference value keeping in mind that it 
may underestimate the number fraction of type Ibc to type 
II supernovae. 
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Z=0.040 >'Z=0.020 
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He mass in ejecta 

Fig. 3. Correlation between the mass of nitrogen in the 
supernova ejecta and the mass of hehum. The SN II region 
(red) and SN lb region (blue) are represented, as well as the 
SN Ic region (on the left of the short dashed line). Masses 
are given in Mq (color figure available online). 



The choice of m'^™ is more delicate. We first note that 
there are no models without any helium in the ejecta. For all 
the metallicities and masses considered here, the minimum 
value of the mass of helium in the ejecta is 0.31 Mq. Second, 
we see that passing from the value 0.4 to 0.5 changes the 
mass limit a lot between the type lb and the type Ic. Still, 
some changes are brought when one passes from 0.5 to 0.6, 
but in general are much more modest except at solar mctal- 
licity. Choosing values for m^^" between 0.6 and 1 - 1.5 M0 
would not change the results significantly when adopting 
0.6. From these considerations it appears that 0.6 is a kind 
of limiting value, because below that value, the mass limit 
between type lb and type Ic presents an high sensitivity to 
the exact value adopted for m^™ (for instance, choosing 0.4 
would make type Ic appearing only in a narrow interval of 
metallicities). Above it the sensitivity of the mass is much 
weaker and the mass limits are much less dependent on the 
exact value chosen. In that respect the value of 0.6 does ap- 
pear to us as the most reasonable choice. Also this choice 
gives a monotonic decrease in the mass limit between type 
lb and type Ic when the metallicity increases and allows 
all WC and WO stars to be type Ic progenitors. Thus, in 
the following we adopt the limiting values corresponding to 
m™" = 0.6. 

We are now able to determine the type of the supernova 
it will give for each model (see column 7 of Table [T]). Table 
m gives for each metallicity the initial mass ranges leading 
to different progenitors, SN events, and stellar remnants. A 
comparison between the mass ranges for the different types 
of progenitors and of SN types shows that all SG stars pro- 
duce a type II SN. The WNL stars generally lead to type 
II SNe except may be those WNL stars on the verge of be- 
coming WNE stars. They present nearly no hydrogen on 



their surface (see the case of 25 M0 model at Z — 0.02). 
WNE stars contain no hydrogen, but enough helium to be- 
come a type lb SN. WC stars end their life as type Ic except 
the less massive ones that contain too much helium to be 
classified as a type Ic SN and thus would explode as a type 
lb. Finally, all the WO stars lead to a type Ic SN event. 

Using Table [H we can build Fig. [2l showing the dif- 
ferent supernova types expected for various initial masses 
and metallicit i es. As for Fig. [Tl we use the models of 
lEkstrom et all (|2008D to complete the figure for Z = 0.0. 
We see that type lb SNe arise from a mass band recov- 
ering the WNL / WNE band (see Fig. [1]) and the lower 
range of WC stars. Interestingly, the estimated initial mass 
of the progenitor of SN 2008ak, which exhibits a transition 
spectrum between type II and type lb, is in the range of 
25 — 30 Mp) for a metallicit y around solar ([Crockett et al.l 
I2OO8I : iPastorello et al.ll2008f) . This agrees with the present 
stellar models. Type Ic SNe cover a much broader range of 
initial masses than type lb. This trend reflects that most 
stars which, at a given stage, in their evolution have peeled 
off their H-rich envelope do not stop at that stage. Their 
evolution drives them beyond up to the stage where most 
of their He-rich layers have been peeled off. 

One of the main uncertainties of our models is the mass 
loss rate. Some recent studies sugge st that this rate cou ld be 
overestimated by a factor of 2 (Va nbeveren et al]|2007( ). We 
briefly discuss here qualitatively how the results presented 
here would be modified if we had considered lower mass 
loss rates. The main differences are : 

— The entrance in the WR phase is delayed, so this phase 
is shorter. 

— At a given metallicity, the minimum initial mass to ob- 
tain a WR star is increased. 

— According to this last point, the minimum initial mass 
providing a type Ibc SN at a given metallicity is also 
increased. 

Globally, the decrease in the mass loss rate in our models 
corresponds to a shift in the limits shown in in Figs. [1] and 
[2]up along the metallicity (vertical) axis. 

6. Nature of the stellar remnants 

Superposed to Figs. [1] and [21 we have indicated, as a light 
grey zone, the regions where a black hole (BH) might be 
formed instead of a neutron star (NS) (see also Table[3|). For 
drawing this zone, we assumed that 2.7 Mq is the maximum 
ma ss of a NS. This value is com patible with the one given 
by IShapiro fc Teukolskvl ([1983f ) and a lso with the recen t 
discovery of a massive NS of 2.1 M© ([Freire et al.ll2008[ ). 
By adopting this mass limit and the relation Mns ve rsus 
Mco deduced from the models of iHirschi et al.l ([2005f ). we 
can estimate the mass ranges of stars producing a NS for 
each metallicity, or a BH. For Z < 0.01, all stars with an 
initial mass higher than 30 Mq finish their life as a BH, 
that is, all WR stars and the most massive red- or blue- 
supergiants. The minimum initial mass to become a BH is 
higher at higher metallicities. There is also a maximum ini- 
tial mass to become a BH at high metallicity. This comes 
from the intense stellar winds that peel off the most mas- 
sive stars and prevent them from having sufficient massive 
cores to produce a BH. As a result, the mass range of BH 
progenitors decreases when more and more metal-rich en- 
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Fig. 4. For each metallicity, the left column represents the 
ratio type lb SNe / type II SNe, and the right one the 
ratio type Ic / type II. The contribution of each kind of 
progenitor is also represented in each column. We used our 
reference case to built this figure (color figure available on- 
line) . 



vironments are considered, and for Z >^ 0.040, no BH is 
expected from the single-star scenario (see Fig. [T|). 

We note that, except at low metallicity {Z < 0.01), the 
majority of WNL and WNE stars produce a NS. WC stars 
may produce either NS or BH for metallicities higher than 
about 0.01. Below that metallicity they would only pro- 
duce BH. WO stars produce only BH. In the same manner 
looking at Fig. ^ we see that most of type Ib's produce a 
NS. Type Ic's produce either a BH or a NS for metallicities 
above 0.02. For lower metallicities, only BH are expected 
to result from a type Ic event. 

7. Frequency of type lb, Ic SNe, and comparison 
with observations 

In Fig. m we show the variation with metallicity in the num- 
ber ratio of type lb, Ic SNe for each type of progenitor with 
respect to type II SNe. To obtain these results, we used a 
Salpeter initial mass function (IMF), and apply the relation 
(in case of the Ib/II ratio): 



SNIb 
SNH 



M„ 



M 



-2.35 



dM 



r-M„ 



(1) 



where Mmin.b is the minimum initial mass giving a type 
lb SN, Mmax. b the maximum initial mass giving a type lb 
SN, and Mmax.ii the maximum initial mass giving a type 
II SN . We used the mass limits given in Table [H and we 
supposed that an SN event occurs, even when a BH forms. 
Similar expressions are used to compute the other ratios 
with respect to type II SNe. 



The first striking point when looking at Fig. [4] is that, 
for most of the metallicities studied here, type Ic are more 
frequent than type lb. This was already apparent in Fig. [2l 
but here it is confirmed after the mass ranges involved have 
been weighted by the IMF. Only at solar metallicities would 
present models predict a similar frequency between lb and 
Ic. A second point deserving to be mention is that the vast 
majority of the type Ibc supernovae come from WC or WO 
stars. WN stars contribute only to a relatively small pro- 
portion of the type Ib's and only at metallicities equal to 
or above solar. 

These frequencies are compared with recent observed ra- 
tios iii_Jj£j5l_and_our models with the observed rat e s give n 
bv ICappellaro et all (I1999D: iPrantzos fc Boissie^ (|2003[) : 
ISmartt et al.l (|2009l ) and iPrietoejaU (2008 15 We also 
add th e predi ction of bina r y star models of Eldridge et al.l 
and of lFrver et all pOOTf ). The gray areas indicate 



the variation of our computed ratio taking the mass limits 
M^y^^ obtained varying m^'" between and 0.5. 

Our models reproduce the observations up to Z = Z(r0- 
At metallicities higher than about 0.02, present models pre- 
dict too few type Ibc supernovae, although in view of the 
big error bars, it is difficult to ascertain that this is a real 
deficit. The increase in the ratio of type Ibc/type II SNe 
with the metallicity naturally arises from the metallicity 
dependence of the stellar winds. At high metallicity stellar 
winds are stronger, making the formation of stars without 
an H-rich envelope that explode as a type lb or Ic supernova 
easier. 

The results obtained including close binary evolution- 
ary scenarios are also able to reproduce the general increase 
in the SN Ibc / SN II ratio with respect to the metallic- 
ity. In these models most of the type Ibc supernovae occur 
as a result of mass transfer in close binary systems. Thus 
we face here the situation where two very different mod- 
els (single stars with rotation/close binary evolution with 
mass transfer) are both able to give a reasonable fit to the 
data. Actually both scenarios probably contribute to the 
observed populations of type Ibc supernovae. However, it 
would be interesting to know their relative importance and 
how their relative importance changes with the metallicity. 
It might be that both scenarios predict different behav- 
iors for the way the frequencies of the type lb and type Ic 
SNe vary as a function of the metallicity. Below we discuss 
the predictions of single-star rotating models and compare 
them with the available observations. 

From the data of iPrieto et al.l (|2008l ). we extract sepa- 
rate observed rates for type lb and type Ic SNe. They are 
represented for the SN lb / SN II rate, and for SN Ic / SN 
II in Fig. [SI We note that, in the high metallicity range, 
the observed rate of type lb SNe is lower than the rate of 
type Ic by about a factor two. The trend deduced from the 
present theoretical models is in good agreement with the 
observed one. 



^ In their recent paper, iPrieto et al.l (|2008l ) analyze a sample 
of 77 core collapse SNe (and 38 type la SNe) in the redshift 
range 0.01 < z < 0.04. This sample is composed of 58 SNe II, 
13 SNe Ic, 3 SNe lb, and 3 SNe Ibc, all with a measurement 
of log (0/H) -I- 12 for th e liost galaxy. 

^ In lMevnet fc Maedeil (|2005l l. we already made this compar- 
ison, although with a less sophisticated method to link the pro- 
genitor to the supernova type it produces. 
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Fig. 5. Rate of SN Ibc / SN II if all models produce an SN 
(solid line) or if models producing a BH do not explode 
in an SN (dashed line). Grey areas are the correspond- 
ing estimated error s from our models. P entagons are obser- 
vatio nal data fromlPrietq et al.1 (|2008D . triangles are data 
from iPrantzos fc Boissieij "( 20031 ). circle, and square, are 
measu reme nts at solar metallic ity from ICappellaro et aP 
(|l999t) and ISmartt et all (|2009t) re spectively. The dotte d 
line represents the binary models of lEldridge et al.l (|2008t) . 
and the dotted -dashed line the r ate obtained with the bi- 
nary models of iFrver et al.l (|2007t ). 



If we adopt mg^" = 0.4 instead of 0.6 (see Fig. [7]), we see 
that our models never produce type Ic SNe (except a small 
number at Z = 0.4 Zq). In that case, models would predict 
that most type Ibc supernovae would be type lb, which is 
contrary to what is observed. This numerical experiment 
illustrates the strong dependence of the present results on 
the maximum quantity of helium in the ejecta that can be 
accommodated in a type Ic supernova event. In case it were 
confirmed that this value is as low as 0.4 Mq, then present 
models would not provide a good fit to the observed data. 
Enhancing this quantity to values superior or equal to 0.55- 
0.6 suffices to provide a good fit. 

How do the results change when one makes the hypoth- 
esis that when a BH is formed no SN event occurs? The sit- 
uation for the variation with metallicity in the number ratio 
of Type Ibc to type II SNe is shown in Fig. \5\ In Fig. [51 we 
show what would be the predictions of the present models 
for the frequencies of type lb and type Ic taken separately 
and adopting m^™ = 0.6. At low metallicity, no type Ibc 
SNe are expected. At Z 0.010 = 0.5 Zq, type lb SNe 
begin to come out of the BH domain, making the SN lb 
/ SN II ratio begins to increase. For Z > 0.02, single-star 
models might still account for more than half of the type 
Ibc supernovae. However, aX Z — 0.02, all the SNe coming 
from single-star evolution would have type lb, which is not 
consistent with the observations. At twice solar metallic- 
ity, all our models lead to a neutron star, and the ratios 



Fig. 6. Rates of SN Ic / SN II (blue long-dashed line), SN 
lb / SN II (red short-dashed fine), and SN Ibc / SN II 
(black solid line) for the refer ence ca se. The points are ex- 
tracted from the data of Prieto et al.l (pOOSV triangles (red) 
represent the SN lb / SN II rate, and upside-down trian- 
gles (blue) the SN Ic / SN II rate. Each triangle corresponds 
to a sample of 11 core collapse SNe (color figure available 
online). 



are identical to the one obtained in Fig. 6. These compar- 
isons show that, considering that the hypothesis "when a 
BH is formed no SN event occurs" makes a big difference 
for metallicities below about 0.02. Above that metallicity, 
the results are only weakly affected. 

If we take a lower maximal mass for neutron stars, this 
situation becomes more extreme: almost all the WR stars 
produce a BH, thus the rate of SN Ibc / SN II is zero or very 
low at every metallicity. This illustrates the high sensitivity 
of the above results on either the mass limit for forming BH 
and/or the possibility or impossibility to have a SN event 
when a BH is formed. 

Probably, the hypothesis according to which no super- 
nova event appears when a BH is formed is too restrictive. 
Until no clear understanding of the physics involved in core 
collapse supernova explosion is reached, it is difficult to 
draw firm conclusions at that point, especially when one 
considers the possible effects of rot ation. In th at context, 
the collapsar scena rio for GRB s (Wooslev 1993) needs the 
formation of a BH (jPessart et al. 2008,) . and this formation 
is at least accompanied in some cases by a type Ic super- 
nova event. Also, the observation of the binary system GRO 
J1655-40 containing a BH (Israelian et al. 1999) suggests 
that a few stellar masses have been ejected when the black 
hole formed and therefore that an SN event occurred. This 
is suggested by the important chemical anomalies observed 
at the surface of the companion whose origin is attributed 
the companion being accreted a part of the SN ejecta. These 
examples show that there are at least some circumstances 
where the formation of a BH does not prevent an SN event 
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Fig. 7. Same as Fig. H but with = 0.4 Mq instead of 

0.6 (see text). 
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Fig. 8. Same as Fig. [51 with the assumption that when 
a black hole is formed, no SN event occurs (color figure 
available online). 



to occur. Probably, reality is somewhere between the two 
extreme cases discussed before: "all massive stars produce 
an SN, independently of the remnant type" , and "only the 
stars leading to a neutron star produce an SN" . 

As in Sect. [HI it is interesting here to infer how a decrease 
in the mass loss rate modifies the previous results. At first 
order, a decrease in the mass loss rates is equivalent to a 
shift in the curves presented in Figs. [5] to [S] to the right 



along the metallicity (horizontal) axis. In general, in that 
case, the theoretical ratios will be lower than the observed 
values. 



8. Gamma ray burst progenitors 

It is now admitted that there is a connection between 
so-called long-soft GRB events and type Ic supernovae 
(IWooslev fc Bloomll2006[ ). As seen in the previous section, 
this kind of supernova can be produced by both WC or 
WO stars. In this way, WC and WO stars are natural can- 
didates to be GRB progenitors. Howe ver, GRB are primar - 
ily found in metal-poor environments (iModiaz et al.|[2008l ). 
while type Ic's SNe appear mainly at high metallicity. Many 
physical reasons have been invoked to explain why GRBs 
seem to occur only at low met allicities. Among them a re 
the following (see the review bv IWooslev fc Bloom|[200l) : 

— At low metallicity, stellar winds (even during the WR 
phases) are weaker, thus bringing away small quantities 
of angular momentum. Black holes are also more easily 
formed since higher final masses are obtained. 

— At low metallicity the transport of angular momentum 
between the core and the envelope is less efficient than at 
high metallicities, because of slower meridional currents 
in metal poor stars. 

— Since the chemical mixing due to rotation is more ef- 
ficient at low Z , homogeneous evolution is more easily 
obtained in metal-poor regions. Homogeneous evolution 
allows massive stars to produce a type Ic SN event with- 
out having to lose large amounts of mass. Indeed, a per- 
fectly homogeneous evolution (actually never realized) 
would allow the formation of a pure CO core (and then 
lead to a type Ic SN event) at the end of the core He- 
burning phase without the need for the star to lose any 
mass! 

— The distribution of initial velocities at low metallicity 
might contain mor e fast rotators than at high metallic- 
ities (see Fig. 9 in iMartavan et al.ll2007[) . 

At the moment, these arguments remain quite speculative. 

What can be done presently is to compare the observed 
GRB frequency with the observed frequency of potential 
candidates. First, as mentioned above, type Ic supernovae 
at low metallicity do appear interesting candidates. From 
Fig. [HI one can see that the observed rate of type Ic super- 
novae from single star models is still above the estimated 
number ratio GRB / core collapse supernovae (CCSNe) 
even when one only considers low-metallicity, type Ic SNe. 
Here we supposed that the formation of a BH does not pre- 
vent an SN event. In case of course only rare circumstances 
would allow an SN to occur when a black hole is formed, 
then the situation might be very different (see below). 

Other interesting candidates are the WO stars that pri- 
marily occur at low metallicity. Table [5] lists all the known 
WO stars. If we take the value log(0/H) + 12 = 8.66 
([Grevesse et al. 2007) for the solar oxygen abundance, we 
see that 6 out of a total of 8 occur in regions with a metallic- 
ity under Z/Zq < 0.9. These stars will explode as a type Ic 
SN. The frequency of type Ic's SNe with WO star progeni- 
tors is shown in Fig. [5] The expected rate is only marginally 
compatible with the observed GRB rate (assuming that the 
aperture angle of the bipolar jet is very small, typically 
around 1°). This conclusion has already been obtained by 
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Fig. 9. Observed rate of type Ic SNe (triangles with error 
bars) and predicted rate of type Ic SNe whose progenitor is 
a WO star (dashed hne) with respect to the total number of 
core collapse SNe (CCSNe). The gray rectangle represents 
the extension in metallicit y (|Modiaz et al.ll2008l ) and rate 
(jPodsiadlowski et al.|[200i ) of GRB events. 



iHirschi et al] (l2005h . so even restraining the progenitors of 
GRB to WO stars would still not match the observed fre- 
quency of GRBs. 

From the above discussion, it can be deduced that ex- 
ploding as a type Ic SN in metal-poor regions (or having 
a WO progenitor) is not a sufhcient condition for a GRB. 
Physical characteristics shared by a subsample of the metal 
poor type Ic events exists that are needed to obtain a GRB 
event. Probably this physical characteristic is the high an- 
gular momentmn in the core (Yoon et al. 2006; Woosley & 
Heger 2006). 

If GRB would only occur for initially very fast-rotating 
stars, rotating so fast t hat the s e sta r s would follow a 
homogeneous evolution (iMaedeil Il987t lYoon et al.l l2006t 
IWooslev k Hegeil 120061: [Mevnet fc Maederl I2007D ." can we 
expect to find any peculiar feature in the chemical com- 
position of the ejecta testifying this previous homogeneous 
evolution? Or in other words, is there any difference in the 
chemical composition of the ejecta between a type Ic having 
"normal WC or WO" progenitors and those arising from a 
model that followed a homogeneous evolution during the 
MS phase? To that purpose we computed the evolution 
of a 60 Mq model followed up to the end of central Si- 
burning, with an initial metallicity of Z = 0.002 and an 
initial rotation rate of n/flcrit = 0.75. T he model was com- 
puted including magnetic field, following ISpruitj ()2002f ) and 
iMaeder fc Mevnet (|2005l ). This model follows a nearly ho- 
mogeneous evolutionary track during the MS phase. The 
angular momentum content at the time of the presuper- 
nova is sufficient for producing a coUapsar. In Tables [1] and 
O we present various characteristics of that model. 



Table [T] shows that the model is in the WC phase at 
the end of its evolution. The mixing induced by the mag- 
netic field produces larger cores than in the models without 
magnetic field. Thus, due to mass loss, the core is uncov- 
ered at an early stage of He-burning, leading to a WC star. 
The helium content of the ejecta (0.64 M©) is also slightly 
greater than in our other models of similar mass and is 
also just above the value of 0.6 Mq below which a type Ic is 
supposed to occur in our classification scheme. This implies 
that our homogenous progenitors lead to a type lb SN. On 
the other hand, as mentioned above, the limit of 0.6 has 
not to be taken as a firm limit, and stars with higher He- 
content might also explode in a type Ic event. Moreover, 
the models having followed a homogeneous evolution will 
end their evolution with more angular momentum in their 
central regions. This may have a significant impact on the 
way the stars explode, as well as on the spectral features 
arising from the SN explosion! 

It is also interesting to compare the chemical composi- 
tion of the ejecta as they are obtained in the present models 
when they follow a homogeneous evolution or not. Compare 
for instance the 60 M© at Z = 0.004 (heterogeneous case) 
with the 60 M© at Z = 0.002 (homogeneous case). The 
differences in the masses of He, CNO elements, and Z are 
very small. Thus there is little chance from the observa- 
tions of the composition of the ejecta in CNO elements to 
be able to distinguish between a normal and a homogeneous 
evolution. 

A point that would be interesting to check in a later 
work is the following: since the star rotate very fast dur- 
ing a great part of its evolution, it will produce anisotropic 
stellar winds (see Fig.[ini), typically bipolar winds and prob- 
ably equatorial mass loss when the critical limit is reached 
(Maedcr 1999). These anisotropic winds will s hape the cir- 
cumstellar environment of the star (Ivan Marie et al.ll2008l 
and see also Fig. [TU]) and it might be that some traces 
of the resulting particular morphology will still be present 
at the time of the SN event (for instance, equatorial mass 
loss probably gives rise to a slow equatorial expanding disk 
whose traces might still be present when the star explodes) . 
In that case the circumstellar environment of GRB may be 
peculiar. This can in turn have an impact on some features 
in the spectrum. More generally, the circumstellar environ- 
ment of stars that have been in a not too remote past very 
fast rotators may also be characterized by such features. 
It would of course be very interesting to find in the cir- 
cumstellar environment of some stars, now slowly rotating, 
traces of their very rapid rotation in a previous phase of 
their evolution. 



9. Conclusion 



In iMevnet fc Maederl (|2005f ) we showed that single-star 
models can account for the increase with the metallicity 
in the ratio (SN Ib-fSN Ic)/SNII. In the present work we 
obtained the following additional results: 

— Single-star models indicate that the mass of helium 
ejec ted in type Ic superno vae is at least 0.3 M© (see 
also lEldridge fc Toutlf200l . 

— Type lb supernovae may be produced by WNL, WNE, 
or less massive WC stars. Type Ic SNe are only produced 
by WC or WO stars. 



12 

Table 5. WO stars catalog. 
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Name Galaxy galactoccntric radius log(0/H) + 12 Cat number type Ref 



DBU 1 


MW 


4.6 


8.76^ 


WR 93b 


W03 


3, 4 


MS 4 


MW 


9.2 


8.56^ 


WR 30a 


W03 


1 


Sand 4 


MW 


3.4 


8.81^ 


WR 102 


W02 


1, 3 


Sand 5 


MW 


7.8 


8.62^ 


WR 142 


W02 


1 


Sand 1 


SMC 




8.05 




W04 


5, 6 


Sand 2 


LMC 




8.35 




W03 


5, 3 


DR 1 


IC 1613 




7.90 




W03 


7, 8 


NGC 1313 #31 


NGC 1313 




8.23 




W03 


9 



(1) va n der Hucht (2001). (2) Rood ct al. (2007) . (3) Drew et al. (2004) . (4) van der Hucht (2006) . 
(5) Hunter et al. (2007). (6) St-Louis ct a^^ (2005') . (7) Kin gsburgh fc Barlow (1995) . 
(8) Bresolin et al. (2007) . (9) Hadfield fc Crowther (2007) 



Note: 1st column is the name, 2nd one is the host galaxy, 3rd one is the galactoccntric radius (only for Galactic WO), 4th one is 
the oxygen abundance, 5th one is the number in The Vllth Catalogue of Galactic Wolf-Rayet Stars^' 6th column is the type of 
the WO star, and the last one the references. 




— In the framework of the "BH-no supernova event" hy- 
pothesis, all type Ibc supernovae at low metallicity 
should come from close binary evolution. At solar and 
highcr-than-solar metallicities, more than half of the 
type Ibc might still be explained by single-star models. 

— Our models predict too high a rate of GRB, even con- 
sidering that only WO stars lead to this rare event. 
Probably GRB only originate in the most rapid rota- 
tors. 

— There is no way to distinguish on the basis of the chemi- 
cal composition of the ejecta a type Ic supernova arising 
from a normal heterogeneous evolution from a type Ic 
supernova arising from a homogeneous evolution. The 
past fast rotation of the GRB progenitor may have left 
some imprint on the circumstellar matter distribution. 

— For Z > 0.04, all massive stars produce neutron stars. 
When the metallicity decreases, the range of initial mass 
producing neutron stars decreases. Below a metallicity 
of about 0.01, all stars more massive than about 30 Mq 
produce a black hole. 



Fig. 10. Young bi-modal shaped nebula blown from the 
highly aspherical, polar dominated winds around a 20 Mq 
star at Z = 10~^. The star has evolved during 21552 yr 
since it has left the ZAMS. The maximal extension of the 
wind is ^ 1.5 • lO^^cm. This model was computed by the 
A-Maze code (jWalder fc Folinill2000f) (color fi gure available 
online) . 



— If we consider that the ejecta of a type Ic supernova may 
contain up to 0.6 Mq of He and that an SN event occurs 
even when a BH is formed, single-star models are able 
to reproduce the observational cumulative ratio of SN 
Ibc / SN II well along with both detailed ratios SN lb / 
SN II and SN Ic / SN II. Because of the small number 
of observed events with known log(0/H), the actual es- 
timations of the ratios SN Ib/SN II and of SN Ic/SN II 
have large uncertainties. Additional measurements will 
permit better constrained stellar models. 



It does appear difficult on the basis of the variation in 
the ratio (SN lb -I- SN Ic) /SN II with the metaUicity to dif- 
ferentiate the relative importance of the single and binary 
channels. As can be seen from Fig. O scenarios with both 
a significant proportion of close binaries and with a mod- 
est one are able to reproduce the observational trend. To 
make progress it would be interesting to study the predic- 
tions of the close binary channel for the relative frequency 
of type lb and type Ic supernovae and to compare them 
with the observations. Also it would be useful to have pre- 
dictions similar to those presented here for what concerns 
the chemical composition of the ejecta. 

One can also note that, in case the hypothesis "no SN 
occurs when a BH forms" is correct, the binary scenario is 
responsible for most, if not all, type Ic and type lb events 
at low metallicities, while at high metallicities single stars 
would still account for more than one half of the type Ibc 
supernovae. In the framework of that hypothesis, the long- 
soft GRB (at least those associated to type Ic events) would 
originate in binary systems. Let us note however that some 
observations show that at least in some circumstances, a 
SN occurs when a BH is formed. 
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A promising way to study the respective importance of 
the single and close binary channel for the formation of type 
Ibc SN progenitors is to look for single-aged stellar popu- 
lations with masses at the turn off between 8 and 25 M0 . 
If such clusters very rarely show both RSG and WR stars, 
then this can be interpreted as these two types of stars 
coming from very different mass ranges just as would be 
expected in single-star scenarios. If on the contrary, the 
simultaneous occurrence of both RSG and WR stars is fre- 
quent, then it may indicate that there is an overlap of the 
initial mass ranges of progenitors of RSG and WR stars. In 
that case, close binary evolution is probably required (ex- 
cept maybe at high metallicity where the overlap also exists 
for single-star scenarios). 
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